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O 

The confining and topological properties of the compact Abelian Higgs model with 
doubly-charged Higgs field in three space-time dimensions are studied. We consider 
i-C ■ the London limit of the model. We show that the monopoles are forming chain- 

ry | like structures (kept together by ANO vortices) the presence of which is essential 

for getting simultaneously permanent confinement of singly-charged particles and 
breaking of the string spanned between doubly-charged particles. In the confinement 
phase the chains are forming percolating clusters while in the deconfinement (Higgs) 
| phase the chains are of finite size. The described picture is in close analogy with 

the synthesis of the Abelian monopole and the center vortex pictures in confining 
■ non-Abelian gauge models. The screening properties of the vacuum are studied by 

means of the photon propagator in the Landau gauge. 

in ■ 

PACS numbers: 11.15.Ha,11.10.Wx,12.38.Gc 

I. INTRODUCTION 

^ ■ The Abelian Higgs model with compact gauge fields (cAHM) has potential applications 

! in various areas of physics. Examples range from high energy physics jH 0] over condensed 
matter physics 0, 0] to neural networks [5] etc. Recently a closer relation of this model 
to the infrared sector of the theory of strong interaction (QCD) was discussed jf|. The 
model is also interesting by itself due to its non-perturbative features which can be treated 
analytically in some cases. 

The non-perturbative features of the cAHM arise due to presence of topological de- 
fects among the excitations of the model. There are two types of defects: monopoles and 
Abrikosov-Nielsen-Olesen (ANO) vortices Q. The monopoles appear due to the compact- 
ness of the gauge vector field while the existence of vortices is linked to the compact- 
ness of the phase of the Higgs field. Vortices carry magnetic flux which begins/ends on 
monopoles/antimonopoles. 

The study of the cAHM with doubly-charged Higgs field (Q = 2 cAHM) presented 
below is motivated by the relation of this model to gluodynamics. In the confinement phase 
of the Q = 2 cAHM the singly-charged external particles must be confined by a linear 
potential while the potential between external doubly-charged particles should flatten at 
large separation between the charges due to the so-called string-breaking phenomenon. This 
feature makes Q = 2 cAHM similar to gluodynamics in which the fundamental charged 
(static quarks) are confined while the potential between the adjoint charges (static "gluons") 
is asymptotically screened. 
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The different effect on singly- and doubly-charged external particles can be reconciled 
with monopole dynamics only if one assumes that monopoles are bound into chains which 
occupy the whole volume of the lattice. Indeed, the monopole gas is known |8[ to be confining 
at all distances, and no string breaking and, consequently, the flattening of the potential 
is possible. As a relevant physical example one may cite the compact U(l) gauge model 
without matter fields (i.e. in compact electrodynamics, CQED3) at zero temperature where 
the monopoles form a Coulomb plasma 0. 

The other extreme, the vacuum filled only by a gas of monopole-anti-monopole bound 
pairs, cannot be confining at all This situation is realized in cQED 3 at high tempera- 
tures [13] as well as in cAHM with singly-charged dynamical matter fields 0, E3| • 

Therefore, in the case of Q = 2 cAHM, the monopoles must form structures different 
both from the monopole gas and the magnetically neutral dipole gas. In fact, it was found 
numerically that in the presence of a doubly-charged Higgs field the monopoles actually 
form chain-like structures. This offers an explanation of both confinement of singly-charged 
electric particles and string breaking for doubly-charged test particles. 

The described pattern of charge confinement in Q = 2 cAHM has a close analogy with 
gluodynamics where tight correlations between Abelian monopoles and center vortices (each 
in the respective Abelian projection 1 ) have been found It also suggests a natural mech- 
anism for the formation of monopole sheets (instead of chains in 3D) in 4D gluodynamics. 
For example, in the pure SU(2) gauge model (chosen here for simplicity of discussion) the 
Abelian monopoles are defined with the help of an Abelian gauge, in which the off-diagonal 
gluons (originally ignored in the Abelian projection) play the role of the doubly-charged 
matter fields coupled minimally to the leading diagonal gluons. The presence of the doubly- 
charged dynamical matter leads to the formation of the vortices carrying magnetic flux |l5[ . 
These "magnetic" vortices are similar to center vortices. The off-diagonal matter fields also 
may cause, as pointed out in Ref. p, the monopole trajectories to be confined inside sheets. 

In this paper we complete the study of the cAHM in three space-time dimensions along 
the direction outlined in Ref. restricting ourselves to the London limit. We hope to 
come back to the model at finite quartic Higgs coupling on another occasion. We study 
the structure of the monopole chains, vortex and monopole cluster properties at various 
couplings of the Q = 2 cAHM model. The (Debye) screening properties of the vacuum are 
investigated with the help of the photon propagator. 

The structure of this paper is as follows. In Section |H] we describe general properties of 
the model on the lattice. In this Section we also discuss the monopole and vortex structures 
and derive their interactions by rewriting the partition function using the Villain form of 
the model in the London limit. The topological nature of the confining properties of the 
model is pointed out. In Section ILTT1 we numerically investigate the phase structure and the 
percolation properties of the monopole chains using the model with the cosine form of the 
action. Both finite and zero temperatures cases are studied. In Section HV1 we numerically 
investigate the photon propagator in Landau gauge crossing the phase transition from the 
confined to the Higgs phase. Our conclusions are discussed in Section IV1 



1 For reviews of the monopole and center vortex confining mechanism as well as for a discussion of the 
corresponding gauges/projections an interested reader may consult Ref. 01 an d Ref. 0], respectively. 
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II. THE MODEL 

A. Lattice formulation 

We consider the three-dimensional Abelian Higgs model with compact gauge fields 9 X)fl G 
[— n, 7r) which are defined on links {x,fi}. Let the scalar Higgs field, Q x G C, be doubly- 
charged and written in the form 

% = PJ^ (1) 
where p x is its modulus (radial part) and <p x its phase. The action of the model is 

S = -(3 ^2 cos 6p ~ K cos (-<^ - 2 9 x>fJ + ip x+fl ) + ^2(p 2 x + X(p 2 x - l) 2 ) , (2) 

P X,fi X 

where 9p is the plaquette angle. The parameter (3 is proportional to the inverse gauge 
coupling squared, j3 = l/(ag 2 ), k is the so-called hopping parameter and A the quartic 
Higgs self coupling. 

In a special limit of the model - the so-called "London limit" - the coupling A is infinite 
and therefore the radial mode of the Higgs field is frozen to unity. In this paper we study 
this limit in detail and relegate the more general case of the radially active model to another 
paper. In the London limit Eq. (J2J) reduces to 

^London [0] = ~P ^ COS 9 P ~ K ^ C0S (-<^ - 2 9 X ^ + if x+i x) . (3) 

P X,fl 

The phase structure of the reduced model can be sketched starting from the following 
considerations [J, Hf. At vanishing hopping parameter, k — 0, the model © reduces to the 
pure compact Abelian gauge theory. On the symmetric lattices this model is confining at 
any coupling (3 due to the presence of the monopole plasma ja|. Thus the Iow-k region of 
the phase diagram corresponds to the confining phase ("confining region"). At large values 
of k (called the "Higgs region") the model reduces to the Z 2 gauge model. Indeed, at 
very large k we get the constraint ^ x +ji — 2 9 x ^ — ip x = 2nl x ^ with l x ^ e 'Z . However, 
in the unitary gauge, all (p x — 0, and the above constraint along with the compactness 
condition for the gauge field, 9 X ^ G [— 7r,7r), gives two possible solutions: 9 X ^ = and 
9 X)l x = TV. One can rewrite the model of Eq. (jHJ) in this gauge as the ZZ 2 gauge model for the 
variables z x li = e i9x ^ 1 = ±1. The ZZ 2 gauge model in three dimensions has a second order 
phase transition at the critical coupling [3f 2 ~ 0.7613 ^6[. A transition line separates the 
confinement phase [closer to origin, (/3, k) = (0, 0)] from the Higgs phase [closer to the point 
(J3, k) = (oo,oo)]. 

At zero value of the coupling constant j3 the model is trivial. At large f3 the model 
reduces to the three dimensional Xy-model due to the constraint for the plaquette variable, 
9p = 2-Knp, with np G /Z . This constraint is solved in the form 9 X ^ = —Xx+fi + Xx + 2ttI x ^, 
where Xx G [— 7r, tt) is a compact scalar field. Then a re-scaling of the spin field Xx by a factor 
of 2 gives us the XV-model which is known to possess a second order phase transition 
at = 0.45420(2) between a symmetric and a Higgs phase. The schematic view of the 
phase diagram is shown in Figure ^ 

From our previous studies in the London limit of the Q = 1 cAHM at (3 = 2.0 we 
know that there is no transition but the crossover from the "confining region" to the "Higgs 
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FIG. 1: The phase diagram of the 3D Q=2 Abelian Higgs model in the London limit. 

region" is signalled by a rapid drop of the monopole density. In addition, a string breaking 
phenomenon affecting fundamental test charges (being in the same representation as the 
fundamentally charged matter fields) has been observed in the Higgs region whereas they 
are bound by a linearly rising (string) potential in the confining region. While the drop 
of the monopole density signals the onset of deconfinement, the anomalous dimension of 
the photon propagator is turning to zero at the same crossover k c (/3). This accompanies 
the transition from the symmetric/confining to the Higgs region, whereas the observed non- 
vanishing photon mass at larger k is simply generated by the Higgs effect. 



B. From Q = 1 to Q = 2 and beyond 

In the compact U(l) model without dynamical matter fields external static charges q 
(measured in units of some electric charge e) are permanently confined in three space-time 
dimensions at zero temperature The mechanism is simple: the compact model possesses 
monopoles which form a Debye plasma state. If the external charges are immersed into 
this plasma moving around a Wilson loop, monopoles form a sandwich-like structure of 
magnetic excess charge of opposite sign on both sides of a certain surface. This surface is 
the minimal area surface spanned by the trajectory of the "electric" charge. The thickness 
of this structure is of the order of the Debye screening length. This leads to a finite excess 
free energy per unit area of the double sheet. Thus, oppositely charged external sources at 
distance R are confined by a linear potential V q (R) due to the monopole plasma. 

The situation becomes more complicated when dynamical matter fields with electric 
charge Q e are present where Q is an integer. In this model - known as the compact 
Abelian Higgs model (cAHM) - monopoles are not the only topological objects present 
in the vacuum. The model additionally possesses so called vortices which carry a certain 
magnetic flux. In the general case the total flux carried by the vortex is $ = 2ir/(Qe). 
According to the Dirac quantization condition the monopole carries the charge m = 2n / e. 
The magnetic flux emitted by an (anti-) monopole is squeezed into the vortices (the Meissner 
effect). The number of vortices attached to a single (anti)monopole is Q = m/$. 



5 



If the matter field carries unit charge, Q = 1, only one vortex is attached to each 
(anti)monopole. When the vortices possess a non-zero "vortex tension" (i.e. action per 
length) the monopoles must form monopole-antimonopole pairs ( "monopoliums" or, in other 
words, magnetically neutral dipoles) confined by a linear potential due to the vortex. The 
monopolium has a characteristic size, i? m , which is inversely proportional to this vortex 
tension. This monopole binding changes now the primordial confining potential V q ^Q(R) 
between static external electric charges q e. 

Firstly, we ignore the Meissner effect and suppose that the vacuum consists of pairs 
of size R m . When the monopolium size is small compared to the distance between the 
external electric charges, R m <C R, the influence of any monopole on the potential V^q(R) 
is almost compensated by that of an antimonopole belonging to the same pair. Therefore 
the influence of any ( ant i-) monopole on the long-distance part of V q< Q(R) is negligibly small 
and confinement is lost. In other words, a vacuum formed out of dipoles is not confining, so 
for R 3> R m the potential V q> Q(R) becomes flat. If R m ^> R, the field of a single monopole 
is practically unscreened by the presence of other monopoles. As a consequence, we may 
expect a piecewise linear potential due to the usual plasma-like mechanism. This picture is 
exactly the string breaking mechanism. While the string breaking has to appear for external 
charges of arbitrary strength q, the typical string breaking size might depend on q. So we 
conclude, that the binding of the monopoles into magnetically neutral pairs may successfully 
describe the string breaking in terms of topological degrees of freedom. 

The so far neglected Meissner effect provides an additional distance scale, which is equal 
to the inverse mass of the compact gauge field, Rg = l/M 9 . This scale defines the width of 
the physical vortex and also the characteristic distance at which the magnetic field of the 
monopole deviates from the Coulomb behavior due to its squeezing into the vortices. Thus 
the presence of the Meissner effect leads to an obvious modification of the condition that 
R is larger than the string breaking distance. Instead of R 3> R m (ignoring the Meissner 
effect) it is sufficient to have R 3> min(i? m , Rq) in order to find string breaking. 

If the charge of the dynamical matter field is bigger than unity, the situation becomes 
more complicated. Consider, for example, the Q = 2 case. It is known [l| that odd-charged 
external test particles must be confined in this case while even-charged ones must show 
string breaking at suitably large distances. What structures of monopoles and vortices 
might correspond to this picture? If the monopoles would form exclusively neutral dipoles 
the string breaking should affect both even- and odd-charged external sources. Thus, the 
monopoles have to form other structures to explain the mentioned behavior. 

Since there are Q vortices attached to each (anti-)monopole, the vortices eventually force 
the monopoles to form magnetically neutral dipole-like states as in the Q = 1 case. If 
Q = 2, in addition monopoles and vortices may form also wire-like structures as depicted 
in Figure |2fa). In the case of Q = 3 the monopoles and vortices may also form various 
structures like "benzene molecules", non-planar nets etc., as seen in Figure I2^b). In the 
Q = 4 case the monopole clusters must be similar to fullerene "cages" (carbon clusters). 

C. Topological interactions 

Let us consider now in more detail the contribution of the topological defects to the Wil- 
son loop, W q> Q, which defines the potential between two oppositely charged static external 
sources, (W q ,o) oc exp{— T V^q(R)}. Here we have assumed that the Wilson loop has a 
rectangular shape, T x R with T 3> R. If no matter fields are coupled to the gauge field, 
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Q = 0, all non-perturbative contributions to the long-range potential would come from the 
monopole fields alone because the monopole fields are long-ranged in this case. In the pres- 
ence of an electrically charged matter field, however, the magnetic fields are squeezed into 
vortices, which carry $ = ^ units of magnetic flux. These vortices will provide the major 
contribution to the potential between external sources at large distances. The contribution 
of a general field configuration to the Wilson loop is W^ q = exp{i q $c}- where $c is 
the Abelian flux going through the contour C. Therefore each closed vortex trajectory jo 
contributes to the Wilson loop a factor 

WC Q = exp{2m^E( Jo ,C)}, (4) 

where K/(jo,C) is the number of vortices piercing a surface spanned by the Wilson loop 
contour C. The integer number iL is an analogue of the linking number between vortex 
trajectories and the external particle trajectory discussed below. 

The contribution of vortices to the Wilson loop quantum average in the form (jlj) may 
asymptotically lead to an area law provided the ratio q/Q is not integer. An example of 
such a behavior is known from SU(2) gluodynamics where the so-called center vortices were 
shown to give a dominant contribution to the confinement of quarks Thus, we expect 
that in the Q > 1 cases the monopoles may form various types of clusters which range from 
simple wire-like structures to more complicated extended three-dimensional structures like 
nets etc. 

The contribution of such vortices to the external particle potential depends on the ratio 
of the external charges, q, to the dynamical charges, Q. Eq. (jlj) states that if ^ e /Z, then 
the contribution is trivial while in the case 4 ^ /Z the vortices actually contribute to the 
Wilson loop average. In the latter the analogy with SU(2) tells us - the Wilson 

loop may show the area law. Thus we conclude 



string breaking , 
no string breaking . 



(5) 
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It is quite instructive to reproduce these expectations by rewriting the partition func- 
tion of the cAHM 3 in the London limit in terms of topological defects using the so-called 
Berezinsky-Kosterlitz-Thouless [l9l Eol ] (BKT) transformation. For the sake of simplicity 
we use the Villain-type form of the action with couplings (3 and k. The cAHM 3 partition 
function in the Villain form has the form 

7T 7T 

Z = I VO I Vip exp|-/3||d^ + 27™|| 2 -K||d^-g^ + 27r/|| 2 |, (6) 

where n and I are integer- valued forms defined, respectively, on plaquettes, C2, and links, 
Ci, of the lattice. Here and below we use the compact notations of differential forms on the 
lattice (for a review see, e.g., Ref. bn brief, the plaquette angle 6 P on the plaquette 

{x,/iz/} is written as (d6) Xfiu = 9 X)fl + 9 X+ ^ U - 6 x+i>}tJ , - & x , v . The gradient of the phase of 
the Higgs field belonging to the link {x,fi} is (d<p) x = <f x +ft — W The Villain couplings 
are related to the original couplings as follows: 

WMslog^M]" 1 , SW =[21og^] _I . (7) 
L h\p)i L ii(«JJ 

Let us apply the BKT transformation [l2|, E^, H^j first with respect to the gauge field in 
Eq. (jUJ). The sum over the integer rank-two form n in © can be rewritten 2 in terms of a 
sum over an integer-valued one-form p and an integer-valued three-form m, 

E = E E ■ (8) 

n(c 2 )£2Z p(cx)eZ m(c 3 )e^ 

using the decomposition 

n = k[m] + dp , where d&[m] = m . (9) 

Here k[m] is an arbitrary integer two- form satisfying the second equation in Q. On the 
dual lattice (i.e., on the lattice which is shifted by a half lattice spacing forward in all three 
directions) the form *m is an integer-valued scalar field (zero-form) constrained (due to the 
second equation in Q) as follows: 

J2*m* x = 0, (10) 

*x 

what is equivalent to dm = or 5*m = (with the co-derivative 5 = *d*). 

The integer-valued form *m represents, if non-vanishing, the presence of monopoles 2 



on the sites of the dual lattice. For example, if exactly one monopole (antimonopole) resides 
in a particular cube c 3 of the original lattice, then the form *m is equal to ±1 at the 
corresponding site *x. Equation (|10|) means that the total magnetic charge on the lattice is 
zero. 

Using the Hodge-de-Rahm identity 1 = <5A _1 d + dA~ 1 5 where A denotes the lattice 
Laplacian, the form k[m] in Eq. (jUJ) can be rewritten as k[m] = d/Sr l 5k[m] + 5/S.~ 1 m. 



2 Here and below we omit volume-dependent constant factors. 
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Therefore the two-form n can be written as n = SA l m + d(A 1 5/c[m] + p) and we obtain 
the gauge invariant plaquette in the form 

dO + 2<im = d6 nc - + 27r5A- 1 m , D - C - = 6 + 27rA _1 6k[m] + 2np , (11) 

where 6 |n,c ' is a non-compact gauge field, i.e. 6f c ' G (— oo,+oo). Moreover, the sum over 
the integer- valued variable p and integration over the compact variable 9 can be represented 
(again up to a volume-dependent factor) as an integration over the non-compact variable 

e n - c -. 

Of course, the BKT transformation is nothing but a change of variables. After the first 
BKT transformation Eq. © reads as follows: 

OO 7T 

Zcx I V9 nc - I V<p exp|-/3||dr x i| 2 -47r 2 /3(*m,A- 1 *m) 

8 *m=0 

-h\\d(p -Q6 n - C - + Q ■ 27iA' 1 5k[m] + 27rZ'| | 2 j> , (12) 

where we have used the basic relations of the differential lattice formalism: d 2 = 5 2 = 0, 
A = d5 + 5d, *d* = 5, and * 2 = 1. In addition, the original one-form I is replaced 
through a shift by I' = I + Qp to make the link-dependent part of the action (proportional 
to k) independent of p. The global condition (fTU|) is implicitly understood in Eq. ()12|) and 
subsequent equations. 

Next, we apply the BKT transformation a second time, with respect to the compact phase 
of the Higgs field. We proceed similarly to the case of the compact gauge field, representing 
the sum over the one-form /' as sums over some two- form j and a zero-form r. The forms j 
and r are related to /' as follows: 

I' = s[j] + dr , where ds[j] = j . (13) 

Again using the Hodge-de-Rahm identity we introduce the non-compact scalar field (p n - c -\ 

dy + 2nl' = dp n - c - + 2n5A- 1 j , <^ nx - = y? + 2nA- 1 5s[j] + 2vr . (14) 

The partition function (|T2"|) is now written as follows: 



OO OO 



Zoc j V9 nc - I Vcp nc - J2 Yl exp|-^||dr x -|| 2 -4 7 r 2 ^(*m,A- 1 *m) 

S *j=0 

-k\\d V n - c - -Q6 n - C - + 2ttA- 1 5(] + Q k[m})\\ 2 } . (15) 



*m(*c 3 )e^ *i(*c 2 )e^ 
5*m=0 5*j=0 



The integer-valued one-form 

*jm = Q*k[m] + *j (16) 

represents, if non-vanishing, the vortex line defined on a dual link, similarly to the monopole 
zero-form *m. The constraint 

5* Jm = Q*m (17) 
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(whereas 5*j = 0) indicates conservation of magnetic flux: vortices carry a fraction 1/Q of 
the total magnetic flux emanating from a monopole or absorbing by an anti-monopole. 

It is self-evident that the definitions, introduced in passing, of monopoles *m and vortices 
* j m are identical to those obtained directly from the gauge invariant plaquette P = d6 + 2im 

dP 

m = = dn , dm = (18) 

2ir 

and the link L = dip — Q 9 + lul 

jm = ^J-QJ- = <}/ + Q n , dj m = Qm. (19) 

Performing the Gaussian integration over the non-compact fields 3 in Eq. (|15|) we rewrite 
the partition function in terms of monopoles and vortices: 

Z= E E exp{-S d (*m, *j m )}. (20) 

<5*m=0 S*j m =Q*m 

The action of the topological defects is 



S d ( *m, * Jm ) = An 2 p( *m, *m) + An 2 k( *j m , *j m ) , (21) 

~ 1/2 

where M = Q(k/(3) is the tree-level mass of the gauge boson. 

Let us now consider the contribution of the vortices and monopoles to the potential 
Vg t Q(R) between a pair of external particles with charges ±q, separated by a distance R. 
The potential is given by the quantum average of the Wilson loop defined as 



7T 7T 

(w q {J)) Q = ^jvejv v E 



X 



exp{-/3||d0 + 2nn\\ 2 - K,\\d<p - Q6 + 2tt/|| 2 + iq(J,9)} . (22) 

The external current satisfies 5J=0. 

Following the transformations which led us from Eq. © to Eq. (|20jh the t/.e.a of the 
Wilson loop can be rewritten as 

(W q ) Q = (W q ) ph (W q ) d . (23) 

Here 

W^expj-iw^^J)} (24) 

is the perturbative self-interaction contribution of the external loop given by exchange of a 
massive photon, while the non-perturbative factor is due to the topological defects: 

(W q ) d =^- Yl E exp{-S d (*m, *j m )-S int (*m, *j m ; J)} (25) 

d *m(*c 3 )e^ *j m (*c 2 )eZS 

5*m=0 S*j m =Q*m 



3 This can be done in the simplest way by gauge-fixing the Higgs field angles such that d<p n ' c ' = 0. 
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with 

S int = -2iri±(5j m , ^j^J) + 2m^L( *j , J) . (26) 

The first term in Si nt consists of a Yukawa-type interaction between the vortices and the 
external charged particle. Since the vortices and the monopoles are related to each other by 
the constraint (jl7j) . the first term in S- m t also represents an indirect interaction between the 
monopoles and the external charged particle. 

The second term is of topological nature. It is given by the linking number 

H*Jo,J) = (5 3o ,A- 1 J)e% (27) 

between the external particle trajectory J and the one- form *jo- The form *jo is derived 
from the vortex ensemble *j m by subtracting Q*k[m], cf. (|TH|) . such that 5*j = 0. In 
contrast to *j m , the modified vortex ensemble *jo consists of oriented and closed loops. The 



linking number between *jo and J corresponds to the Aharonov-Bohm (AB) effect (21 . 

In the cQED-like limit, M — > 0, the interaction term (|26J) reduces to the usual cQED-like 
interaction between monopoles: 

lim S int = -2iri^-{5j m , A" 1 J) + 2m^-E( *j , J) 

= 2m^(5(j - j m ), A' 1 J) = -2m(k[m], A" 1 J) , (28) 

where we have used Eq. ()27j) and Eq. f)16j) . The interaction ()28|) does not depend on the 
shape of the Dirac line, */c[m], and does depend on the monopole magnetic charges and 
mo nopole locations, *m = 5*k[m\. 

Note that the separation of a general vortex ensemble, *j m = *jo + *j' m , into closed 
vortices, 5 *j = 0, and open ones, S *j' m = *m, is ambiguous. In the sum (|25p the ambiguity 
disappears since all ambiguous separations enter Eq. ()25j) with the same weight. 

The contribution of vortices to the inter-particle potential V q ,Q(R) is twofold. First, the 
closed vortices *jo interact with the electrically charged particles via the AB effect if the 
ratio q/Q is not integer. Secondly, the vortices influence the monopole dynamics via the 
constraint (fTTj) . contributing in this way to the potential indirectly, via the monopoles and 
their interaction with the test particles. 

Actually, orientable vortices with closed magnetic flux lines cannot lead to confinement 
of electric charges 0]. In that paper numerical studies for the non-compact version of the 
AHM3 have shown that confinement is not realized. This is a model in which orientable 
vortices are present (and consequently, an AB effect might be at work) whereas magnetic 
monopoles are absent. On the other hand, a long-range interaction between vortices and the 
electric Wilson loop may be provided only by the AB effect in a form visible in Eq. (|25|). In 
the model presently under consideration the linking of the non- orientable vortices with the 
external electric current is formally expressed by (J27|) in terms of the redefined closed and 
orientable loops *jo- Making this point, we are implicitly leaning upon that we are working 
in a model with a mass gap in which any interactions - except for topological ones like the 
one mediated by the linking number - must be suppressed. 

In four dimensional QCD the so-called center vortices are known to lead to the area law 
of the Wilson loop via the AB interaction [18]. They are providing a dominant contribution 
to the confinement of quarks if they are localized by means the so-called Maximal Center 
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projection. There is no contradiction between the roles of the AB effect in QCD and the 
non-compact AHM since in the QCD case the center vortices are non-orientable worldsheets 
whereas in the case of the 3D non-cAHM they are orientable lines. Moreover, considered in an 
Abelian sense, the non-orientable center vortices can be represented as piecewise orientable 
vortices with monopole trajectories separating segments with different flux orientations from 
each other. The required non-orient ability implies the existence of monopoles living on the 
vortex sheets. 

We will see that all percolating loops in the case of the model studied here (which are 
responsible for confinement) really form an alternating monopole-antimonopole chain which 
makes them physically non-orientable. Thus, we expect that in cAHM 3 the vortices play 
a role to accomplish both confinement and string breaking phenomena, for the respective 
sources, by restructuring the monopoles configurations. 

III. PHASE STRUCTURE IN THE LONDON LIMIT AND THE 
PERCOLATION PROPERTIES OF MONOPOLE CHAINS 

A. The finite-temperature case 

In this subsection we report on our numerical study of the confining properties as seen 
by static external particles with different electric charges q in the case of the CAHM3 with 
doubly-charged matter. In the confining phase of this model only external particles with 
charge q = 1 must be asymptotically confined while particles with charge q = 2 should 
experience the flattening of the potential at large distances. 

To observe this effect we have simulated the model first on lattices 32 2 x 8. The choice 
of the asymmetric lattice (finite T) was not a matter of principle but simply dictated by the 
fact that in the case of symmetric (L 3 ) lattices the potential can be measured only using 
Wilson loops which are notoriously unsuitable for an observation of the string breaking 
effect. We have performed simulations at fixed gauge coupling constant, (3 = 1.2, the 
choice of which was motivated by visualization reasons. Indeed, in order to clearly see the 
monopole structures, the density of the monopoles must be neither too high nor too low. 
According to Eq. (}2"Tj) a singly-charged (anti) monopole enters the partition function with 
the action Sm = 47r 2 ( 5A~ 1 (0) (with A _1 (0) ~ 0.25 in the thermodynamic limit). Thus an 
increase (decrease) of the gauge coupling (3 constant suppresses (enhances) the number of 
monopoles. The chosen value of (3 is ideal from the point of view of our purposes. We have 
used about 10 5 configurations for each value of the hopping parameter n. Some results on 
the finite-temperature case (partially included in this subsection) were previously reported 
in Ref. 0. 

To probe the phase structure of the model we first study the vacuum expectation values 
(v.e.v.'s) of the q = 1 and q = 2 Polyakov loops, 

L g (x) =exp{^^0 4 (£,*)}. (29) 
t=i 

Their expectation values are shown in Figure Efa) as a function of the hopping parameter k. 
At small (large) k the f.e.v.'s of both loops are low (high) corresponding to the confinement 
(Higgs) phase. One can clearly see that the (outside the Higgs phase) confined q = 1 external 
charges more sensitively react to the onset of the Higgs phase compared to the q = 2 charges. 
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FIG. 3: The q = 1,2 Polyakov loops (a) and their susceptibilities vs. k at /3 = 1.2 on a 32 2 x 8 
lattice. The insert in Figure (b) shows the behavior of the q = 1 susceptibility in the vicinity of 
the phase transition with the fit according to Eq. ()3U|) . 



This is a typical signature of the q = 2 string breaking. The abundantly available Q = 2 
particles screen q = 2 external charges which leads to the low sensitivity of the q = 2 loops 
to the phase transition. 

In Figure Efb) we show the Polyakov loop susceptibilities vs. k. The q = 2 susceptibility 
does not indicate any significant change at the transition point while the q = 1 susceptibility 
shows a sharp peak signalling the presence of a phase transition. We have fitted the q = 1 
susceptibility in the vicinity of the maximum by the function [ill ]: 

^= (r4 _ ( Cl (30) 
(C7 2 + (k - « c ) ) 

with 6*1,2, 7x an d k c being the fit parameters. From the fit (shown in the insert of Figure Efb) 
as a solid line) we obtain that the transition happens at 

= 1.300(1). (31) 



To fit the susceptibility in the vicinity of the critical value we have chosen the region such 
that x 2 /d.o.f. ~ 1. This choice applies to all fits of susceptibilities performed in this paper. 

The static potentials between g-charged external particles can be expressed via the 
Polyakov loop correlators, V q! Q(R) = —L^ 1 \og(L q {0)L*(R)) . As it was shown in Ref. 0, in 
the confinement phase the potential for the q = 1 external charges is linearly rising at large 
distances while q = 2 potential shows a rapid flattening corresponding to dynamical parti- 
cles popping up from the vacuum to cause string breaking. In the Higgs phase all potentials 
suffer flattening due to deconfining nature of this phase. The corresponding potentials are 
shown in Figure 0] This result completely agrees with the general form of the phase diagram 
shown in Figure H 

As we have mentioned above, the cAHM contains two types of topological objects - 
monopoles and vortices. The simplest characteristic of a topological defect is its density. 
The monopole and the vortex densities are (the volume is V = L 2 L t for finite temperature 
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FIG. 4: The q = 1,2 potentials in the confinement (k = 1.275) and the Higgs (re = 1.325) phases 
at (3 = 1.2 (from Ref. 0). 



and V = L 3 at T = 0) 

Pmon = y ^ ] | m \ ; Pvort = ^ ] I Jm| j (32) 

*C3 *C2 

respectively. The monopole charge is defined in the standard way (see also (|18)l). 

m = ^d[d6] 2w , (33) 

Z7T 

where [• ■ -] 27r /(27r) denotes the integer part modulo 2tc. Following Ref. [22], the vortex 
current is defined as (cf. (fl^|l) 

j = J-(d[d^+ ^-g[d^] 2 j. (34) 

According to Figure EJa) the densities (exposed in lattice units) of both monopoles and 
vortices are gradually decreasing functions of the coupling k what is partly a finite temper- 
ature effect. In the confinement (Higgs) phase the density of both vortices and monopoles is 
high (low), in agreement with the expectation that the confining properties of the cAHM 3 
are due to topological defects. Both susceptibilities are peaked at the phase transition point. 
The fits by Eq. ()30|) of the monopole and vortex susceptibilities give the following transition 
values of the hopping parameter, respectively: 

K ™ n = 1.301(3) , /t™ rt = 1.302(2) . (35) 

These values are consistent with each other as well as with the critical value obtained with 
the help of the q — 1 Polyakov loop susceptibility discussed above. 

Obviously, in the confining phase the monopoles cannot form a monopole or dipole plasma 
because in these cases both external charges q = 1,2 would be simultaneously confined or 
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FIG. 5: (a) The monopole and vortex densities (from Ref. 6]) and (b) their susceptibilities, vs. k 
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the vicinity of the phase transition. The solid lines show the fits by Eq. ()30j) . 



deconfined. Thus the only possible kind of monopole configurations - which can explain 
both the linearly rising potential for the q = 1 electric charges and the string breaking for 
the q = 2 charges - is a monopole chain schematically plotted in Figure Efa). In a monopole 
chain the monopoles and antimonopoles are mutually alternating. Thus the magnetic flux 
coming from a monopole inside the chain is separated into two parts, squeezed into the 
vortices and, as a consequence, forming a non-orientable closed magnetic flux. Each piece 
of such a vortex carries (in average) a half flux, $ = ^ = it. If such a flux pierces the 
g-charged Wilson loop it provides a (multiplicative) contribution to the loop close to (— l) 9 . 
This leads to a disorder for odd charged external particles (eventually leading to a linearly 
rising potential) whereas particles with even charge are not confined. The described picture 
is very similar to what is observed jl^ in QCD for the correlations between the Abelian 
monopoles and the center vortices. 

In Figures Efa-d) we show typical monopole configurations observed in our numerical 
simulations. In the confining phase the monopoles form chains which either wrap around 
the temperature direction (Figure a)) or fill all the lattice volume (Figure Efb)). In the 
last case the monopole chains are said to be percolating. In the Higgs phase the monopole 
chains wrap only around the temperature direction (Figure Efc,d)). In the next subsection 
we investigate the percolation properties of the monopole chains in the zero-temperature 
model. We want to demonstrate that the mechanism causing the connning/deconfining 
effects is the same as for the finite temperature lattice dealt with in this subsection. 



B. The zero-temperature case at fixed gauge coupling (3 = 2.0 

We report now the study of the doubly-charged cAHM at zero-temperature, on the lattice 
16 3 . From 2 • 10 3 to 8 • 10 4 measurements for every /3, k pair have been collected. First we 
map out the phase structure along a vertical line at (3 = 2.0 in the phase diagram Fig. ^ 
varying the hopping parameter k. We mainly concentrate on the quantitative analysis of 
the monopole chains and their components, vortices and monopoles. 
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(c) (d) 

FIG. 6: Examples of the monopole configurations in the confinement phase at k = 1.275 (upper 
row) and in the deconfinement (Higgs) phase at k = 1.325 (lower row). The (anti-)monopoles are 
shown by the dots while the vortices are shown by the solid lines. 

Similarly to the finite-temperature case, for orientational purposes, we measure the ex- 
pectation values of the Polyakov loops which also in the zero-temperature model is defined 
as a Wilson line closed via the boundary conditions. The expectation values of the q — 1, 2 
Polyakov loops are shown in Figure 13(a). As in the finite T case the v.e.v. of the q = 1 
Polyakov loop is noticeably larger in the Higgs phase (larger k's) and almost constant in 
the confinement phase (smaller /t's). The q = 2 Polyakov loop is almost constant crossing 
the transition. Consequently, the q = 1 susceptibility of the Polyakov loop is peaked at 
the transition point for moderate lattice volume while the susceptibility of the q = 2 line 
is insensitive to the transition. The fit of the q = 1 susceptibility by Eq. (J3*U|) is shown in 
the insert of Figure E^b) as a solid line. The corresponding critical value of the hopping 
parameter is 

k c = K q c =1 = 0.687(1) . (36) 

To clarify the role of the topological defects we plot the density of defects in Figure |Sf a). 
One can clearly see that in this region of the phase diagram the density of monopoles is 
much smaller than the density of vortices. This is not unexpected since at (3 — 2.0 the 
action of an isolated monopole is Sm = 47r 2 /3A~ 1 (0) ps 20, therefore the monopoles must 
be suppressed drastically. One may therefore conjecture that at (3 = 2.0 the leading role at 
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the transition should be played by vortices and not by the monopoles, in other words, the 
dynamics of the monopoles is driven by the vortex networks. 

In Figure E^b) we show the susceptibilities of the monopole and vortex densities. The 
vortices are sensitive to the phase transition while the susceptibility of the monopole density 
does not show any noticeable peak at the critical value of the hopping parameter. The fit 
of the susceptibility of the vortex density by function (J3U|) gives the pseudo-critical value of 
the hopping parameter 

K ™ rt = 0.684(3) , (37) 

consistent with the value obtained from the susceptibility of the Polyakov loop. 

To characterize quantitatively the behavior of the clusters we plot in Figure Efa) the 
distribution Dy of the lengths L of the vortices inside the (mutually disconnected) clusters. 
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The distributions for various k are normalized to unity. To see the difference between the 
phases, we plot the distributions for k values in the vicinity of the phase transition. One 
can see that in the confinement phase the vacuum is populated by two different types of 
clusters since the distributions have two peaks. The vortex lengths in the first type of 
clusters is typically short. The vortices of this type are called ultraviolet vortices. The term 
"ultraviolet" is borrowed from 4D QCD which also possesses Abelian monopole trajectories 
(defined in an Abelian projection) which fall into similar classes of loops. In QCD the 
ultraviolet monopoles can be described by a formula DVY ~ L~ x where A is a positive 
number which does not depend on the lattice volume 24]. This experience leads us to 
conjecture a similar behavior in the 3D cAHM as well, except that the loops are now formed 
by vortices. 

The other population of vortices consists of the so called "infrared" vortices which form 
the second peak in the distribution Dy. The average number of the length L of the vortices in 
this peak is proportional to the number of links on the lattice, 3V, multiplied by an "infrared 
density" pj R which may be also referred to as the "vortex condensate". The vortices from 
this condensate occupy uniformly all the volume of the system and their typical length 
is equal to pY R 3V. The fluctuations of the length of the infrared clusters are expected 
to be random. Therefore, they may be described by a Gaussian ansatz (here we refer to 
Ref. [25| where a similar observation was made for monopole trajectories). It is clear that 
the vortex condensate exists only in the confinement phase (signalled by the infrared peak 
of the histogram) while it is absent in the Higgs (deconfinement) phase 

A similar behavior is also observed for the distributions Dm with respect to the monopole 
number N inside the mutually disconnected vortex clusters. Monopoles are pointlike and 
connected by vortices 4 . These distributions are plotted in Figure EJb). One can also see 
the two-peak structure in the confinement phase (the infrared and the ultraviolet monopoles 
corresponding to long and short vortex clusters) and the disappearance of the infrared peak 
in the Higgs (deconfinement) phase. 



4 Inspecting the vortex clusters without monopole-antimonopole pairs in the confinement phase we notice 
that they belong to the ultraviolet component (loops of short length) which is irrelevant for topological 
(AB) confinement. 
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From Figures EJa,b) we conclude that in the confinement phase a vortex condensate is 
formed. The percolating vortices are carrying the monopoles which are responsible for the 
confining properties of the vacuum. 

In order to characterize the cluster distributions of vortices and monopoles D VjM we fitted 
the distribution functions by 

Dy(L) = C x ,v L~ Xv + C 2y exp - p v IR ■ 3V) 2 } (38) 

and by an analogous form for the monopole distributions 

Dfj(N) = C hM N- x ™ + C 2 , M exp{-a M (iV - pf R ■ V) 2 } (39) 

to determine the condensates p^ R l ' . Equations ()38|) and (|3*9*j) describe the distribution of the 
topological defects as the sum of an ultraviolet component (first, inverse power term) and 
an infrared component (second, Gaussian term). 

As we have seen, the second (infrared) peaks in the monopole number and vortex length 
distributions are absent in the deconfinement phase. Therefore we fit the distributions by 
Eqs. (I38|) and (J39j) only in the confinement phase. The fits of the vortex length distributions 
are shown in Figure ITDT a). One can see that the fitting function describes well the data both 




FIG. 10: Fits of (a) the monopole and (b) the vortex histograms by Eqs. (|HH|) and ftffljl . 

in the low-L and high-L regions. To fit the monopole distributions we have set C^m = since 
there is not enough data to fix the ultraviolet part of the distribution. We have performed 
the fit for iV ^ 6 as it is shown in Figure ITuT bb 

The densities of vortices and monopoles in infrared clusters are shown in Figure ^2 as a 
function of k in the confining phase. One observes that the density of the infrared part of 
both types of topological defects drops down when the phase transition point is approached, 
in agreement with our expectations. Deeper in the confining phase the infrared densities are 
practically saturating the total densities. 

We have fitted the infrared densities of both types of topological defects by functions of 
the form 

pVf(K) = B VM ■ (k v c > m - k)i V ' M , K < K C , (40) 
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FIG. 11: The densities of topological defects in infrared clusters as given by fits (|38jl and (|39|) . The 
fits by Eq. (|40jl are shown by dashed lines. 



assuming a "critical behavior" near the phase transition. Both fits are shown in Figure (fTTJl 
as dashed lines. The transition points are: 

«f = 0.680(3) , k v c = 0.684(1) (41) 

from the infrared monopole and vortex densities, respectively. They are in agreement with 
each other as well as with the values determined by from the susceptibilities of the Polyakov 
loops and the vortex density. The "critical exponents" are different: j M = 0.28(2) and 
7 y = 0.37(1). 

Yet another interesting aspect of the cluster structure is the cluster multiplicity dis- 
tribution D(N c i) which represents the total number of vortex clusters per configuration 
(irrespective of the cluster size). We show D(N c i) for various values of hopping parameter 
k in Figure IT21 The distribution of the cluster number can be well described by a Gaussian 
function, 

D(N cl ) = Ae Xp { J N "-f" )2 }. (42) 

The fits are shown in Figure ^] as solid lines. 

The fitted average cluster multiplicity N d and the width of the distribution A are pre- 
sented in Figure IT31 as functions of k. Both parameters are peaked near the phase transition 
indicating that at the critical value of the hopping parameter the few large infrared clusters 
- existing in the confinement phase - decompose into many smaller clusters as the system 
goes over into the deconfinement phase. According to Figure Ufa) the monopole and vortex 
densities drop slightly with increasing k at k c . Therefore, the increase of the cluster num- 
ber is associated predominantly with a cluster-restructuring process and not with a cluster 
creation phenomenon. Moreover, the increase of the width A indicates the existence of large 
fluctuations in the vortex system at k ~ k c . 



C. The zero-temperature case at fixed hopping parameter k = 1.0 



In this subsection we map out the phase structure along an horizontal line through the 
phase diagram of Fig. ^ At k as large as k — 1.0 the confinement phase is realized only for 
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FIG. 12: The multiplicity distributions of vortex clusters and the corresponding fits by the Gaussian 
distribution Q42|) for various k values. 
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relatively small (3. In Fig. El we show the v.e.v.'s of the q — 1 and q = 2 Polyakov loops 
vs. (3. Only the susceptibility of the q = 1 Polyakov loop shows a peak which leads to an 
estimate 

P e = Pr 1 = 1-410(2) (43) 

using a fitting function x &t (P) analogously to (J3UJ). The position of the phase transition 
defined from the fit of the q = 1 Polyakov loop is shown by a vertical dashed line. 

In the confining phase we have now a density of monopoles comparable in magnitude to 
the density of vortices. In Fig. we show both densities and the respective susceptibilities. 
Values as f3 = 2.0 which have led to a vortex-driven transition described in the previous 
subsection are now belonging to the deconfinement phase. The peak of the susceptibility of 
the vortex density is still more pronounced at the transition, but now the susceptibility of 
the monopole density also shows a little peak. 

The positions of the peaks for the monopole and vortex susceptibilities are slightly shifted 
compared to the peak of the q — 1 Polyakov loop susceptibility, Fig. H31 The fits of the 
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FIG. 15: (a) The densities of the topological defects and (b) their susceptibilities vs. (3. The insert 
in Figure (b) shows the behavior of the monopole and vortex susceptibilities in the vicinity of the 
phase transition with their fits and critical /3's. 



monopole and vortex susceptibilities in the vicinity of the phase transition give the following 
values for f3 c corresponding to the transition, respectively: 

/3 c mon = 1.398(1) , /3 c vort = 1.401(1) . (44) 

One can see that the couplings f3 c obtained with the help of the monopole and vortex 
susceptibilities are very similar to each other. 

Histograms of the vortex length an monopole number distribution for the individual 
vortex clusters are shown in Fig. ^] for various /3 values near (3 C . Fits corresponding to the 
two-component forms (J3~%|) and (J3~9"|) are presented in Fig. [T7| 

The fitted values of the infrared monopole and vortex densities are shown in Fig. ITHI 
together with the total densities as functions of (3. Again deeper in the confining phase the 



22 



— i — i — i i i 1 1 1 1 1 — i — i i i 1 1 1 1 1 — i — i i i 1 1 1 1 

monopole k=1.0 




: 

if i l : 



- 1 — i — i i i 1 1 1 1 1 — i — i i i 1 1 1| 1 — i — i i 1 1 1 1 1 1 — i — i ■ ■ ■ i ■ | 1 — i — i— r 

vortex K=1.0 



000 




10 100 1000 i 

(b) 



FIG. 16: Histograms of the (a) monopole and (b) vortex distributions. 
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FIG. 17: Fits of (a) the monopole and (b) the vortex histograms by Eqs. (|35|) and 1)39(1 . 
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infrared densities practically saturate the total densities. The respective /3 c -values mark- 
ing the transition from the infrared clusters are found fitting the /3-dependence with an 
ansatz pf^i^) similar to Eq. (|40p. The fit curves are shown in the Figure as well. The 
corresponding critical values extracted from the infrared monopole and the vortex densities 
are: 

= 1-424(3) , /3 c y = 1.399(2) . (45) 

The "critical exponents" are rather different again: r y M = 0.86(1) and ^ v = 0.499(4). The 
deviation between the values of the transition points /3 C obtained from fitting the two infrared 
densities (|35)l and from the f3 c values obtained from the defect density susceptibilities 
is a finite volume effect, caused presumably by an enhanced sensitivity of these densities to 
the volume of the system. 



IV. THE LANDAU GAUGE PHOTON PROPAGATOR AT THE PHASE 

TRANSITION 

Here we want to clarify how the behavior of the photon propagator reflects the properties 
of the two phases. We restrict ourselves to the case of zero temperature (symmetric lattice). 
In order to calculate the (gauge dependent) photon propagator directly, a special gauge has 
to be chosen. In our case this is the minimal Landau gauge. This gauge is defined by finding 
the global maximum of the gauge functional 



links 



£cos(0^) ->max (46) 



x,ii 



with respect to gauge transformations G. Details of the gauge fixing procedure can be found 
in Ref. Q. 

The photon propagator in lattice momentum space is the ensemble average over gauge- 
fixed configurations of the following bilinear in the Fourier transformed lattice gauge poten- 
tial A (with Puiku) = (2/a) sin(7r£; M /L M )) 

D,M = <Ai*A-0 ■ ( 47 ) 
constructed out of the gauge fixed links: 

h, = (^) 1/2 $>p( 2 " E U — } ) *„ = 0,±l,...,±^. (48) 

The lattice equivalent of the same continuum p 2 can be realized by different vectors k which 
eventually would reveal a breaking of rotational invariance on the lattice. 

We consider the general tensor structure with two form factors D and F for the propa- 
gator: 

^ (p1 = ( v _m,) bM + m,£M. (49) 

When the Landau gauge is exactly fulfilled, F(p 2 ) = 0. On the lattice, this is actually 
the case as soon as one of the local maxima of the gauge functional (j4^j) is reached. Since 
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D^vip) is only approximately rotationally invariant, the form factor D(p 2 ) scatters instead 
of forming a smooth function of p 2 . 

We have measured the propagator on a lattice of size 40 3 at (3 = 2.0 as function of the 
hopping parameter k. O(50) independent configurations have been considered outside the 
phase transition region and 0(200) (roughly) independent configurations near the second 
order transition. The autocorrelation times have been estimated by looking both at the 
monopole and vortex densities, p mon and p VO rt- In order to fight with huge autocorrelations, 
we separated the used configurations by 10000 Monte Carlo sweeps. 

As in the case of the Q = 1 AHM we parameterize the function D(p 2 ) by 



where Z, a, m and C are fitting parameters. Their meaning is as follows: Z is the renormal- 
ization of the photon wavefunction, a is the anomalous dimension, m is a mass parameter. 
The parameter C corresponds to a 5-like interaction in the coordinate space and, conse- 
quently, is irrelevant for long-range physics. From previous studies we know that this fit 
ansatz well describes the data in pure confining and Higgs phases separated by a crossover 
or a first order phase transition. 

Treating now the same ansatz, we present in Figs. EH the obtained fit parameters as 
function of the hopping parameter k. The optimal fit parameters have been found by 
minimizing x 2 /d.o.f . Outside the critical region (k c pa 0.68 - 0.70) the propagator data are 
perfectly described by the simple fit ansatz (|5U|l. The quality of the fits is the same as in 
the Q — 1 AHM j26[, therefore we do not show a comparison of the measured data with the 
fits ansatz. 

The propagator fit leads to a minimum of the photon mass (Fig. 119( a)). On the Higgs side 
of the transition the photon acquires its mass via the Higgs mechanism. The noted minimum 
could be explained similar to the Q — 1 model via the interference between perturbative 
mass and Debye mass effects. 

As shown in Fig. lTW b) the anomalous dimension a drops towards k c and becomes negative 
near criticality within the used ansatz. It tends to zero from below at larger k, deeper in 
the Higgs region. The strong wave function renormalization effects in the confining phase 
(Z is much larger than one, see Fig. ITffl c)) weaken towards the transition and Z approaches 
a roughly constant value near one, deeper in the Higgs phase. For completeness, the contact 
term C as function of k is shown in Fig. IT9T d) as well. 

In the critical region of the second order phase transition we observed that the propagator 
data for the smallest available lattice momenta seem to indicate a non-vanishing p 2 D(p 2 ) for 
p 2 — > 0, contrary to the used fit ansatz (as long as the fitted mass remains finite). For the 
worst observed case at k = 0.695 the smallest measured value for p 2 D(p 2 ) lies roughly two 
times above the fit curve (at p 2 ~ 0.025). The noticed deviation between data and fit curve 
is exclusively related only to the smallest lattice momenta. Nevertheless, the obtained values 
for x 2 / d.o.f. do not exceed 1.17 and do not differ significantly from the x 2 values outside 
the transition. More sophisticated parameterizations for the propagator and eventually 
measurements for still smaller momenta (larger lattices) are needed to clarify this problem 
what also might influence the actual values of our fit parameters. This, however, was beyond 
the scope of our present investigation. 




(50) 
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FIG. 19: The fitting parameters for the propagator measured on 40 3 at (3 = 2.0 vs. k. 



V. CONCLUSIONS 

Summarizing, we have found that the compact Abelian Higgs model with doubly-charged 
matter field has a rich structure in terms of topological defects: monopoles and vortices. 
This structure has been outlined first by analytic considerations in the London limit and 
in the Villain representation. In particular, the topological character of confinement for 
mod(q, Q) ^ via the Aharonov-Bohm (linking number) interaction has been pointed out. 

Extending first results, obtained for a suitable finite-temperature regime of the model and 
published already in 0, we could show by numerical simulation in the London limit, that 
the main property of the model, the existence of vortex clusters with monopoles attached to 
them is a generic feature of the model. This structure can explain that in the confining corner 
of the (3-k phase diagram only odd-g external charges are confined while even-g charges suffer 
string breaking. 

In the confinement phase we found that the vortex clusters are percolating. Both the 
length distribution of the vortex clusters and the monopole number distribution kept inside 
the clusters give rise to the definition of two order parameters, the "infrared densities" of 
monopoles and vortices. The infrared (percolating) component of the vortex network and 
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the corresponding order parameters vanish in the Higgs phase. 

For not too large lattice size at T = the phase transition is signalled also by the original 
gauge field variables, using for example the susceptibilities of the q = 1 Polyakov loop on 
finite lattices (directly related to the loss of confinement for charge q = 1 external charges). 
More quantitative description one can get from the susceptibilities of the monopole and 
vortex densities, and finally in terms of the two infrared order parameters. Different parts of 
the confining phase (high f3 and small k on one hand and small (3 and high k on the other) 
may differ in the detailed mechanism, being vortex driven (at low k) or monopole driven (at 
low /3). 

We also have shown that the gauge boson propagator changes its character at the phase 
transition to the Higgs phase where it changes from confining to Yukawa type. There are 
remarkable differences from the Q = 1 AHM. Using our simple parametrization ()5U|) for 
the form factor D(p 2 ) we found that in the Q = 2 case both the critical dimension a, the 
photon renormalization function Z and the contact term C go through minima at the phase 
transition point, which is not the case for the Q — 1 case. Moreover, in a vicinity of the phase 
transition the renormalization parameter Z gets smaller than unity and the propagator can 
be described by a negative anomalous dimension contrary to the Q = 1 case. We were 
not able to see the fitted gauge boson mass vanishing at pseudocriticality (eventually to be 
expected from the second order nature of the phase transition). So, this interesting behavior 
can be confirmed only using larger lattice sizes to access still lower lattice momenta and more 
sophisticated parameterizations for the propagator. 

This study deepens our understanding of the phase transition in this model. It also 
further develops the paradigm for the cooperative role that monopoles and vortices play in 
non-Abelian gluodynamics. The physical picture - observed in the present study - has a 
close analogy with QCD (in 4D) where tight correlations between Abelian monopoles and 
center vortices have already been observed. The non-diagonal gluons are usually ignored 
after the Abelian projection (from the maximally Abelian gauge). In fact, they are charge- 
two matter field with respect to the remaining Abelian gauge symmetry. If these matter 
fields would be retained as dynamical, this might naturally lead to the formation of monopole 
sheets on vortex sheets, which in turn should be responsible for the permanent (in the pure 
gauge model) confinement of the fundamental charges (quarks) and for the flattening of the 
potential between the adjoint charges ("static gluons"). 
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